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We study the counting statistics of ultracold bosonic atoms that are released from an optical 
lattice. We show that the counting probability distribution of the atoms collected at a detector 
located far away from the optical lattice can be used as a method to infer the properties of the 
initially trapped states. We consider initial superfluid and insulating states with different occupation 
patterns. We analyze how the correlations between the initially trapped modes that develop during 
the expansion in the gravitational field are reflected in the counting distribution. We find that 
for detectors that are large compared to the size of the expanded wave function, the long-range 
correlations of the initial states can be distinguished by observing the counting statistics. We 
consider counting at one detector, as well as the joint probability distribution of counting particles 
at two detectors. We show that using detectors that are small compared to the size of the expanded 
wave function, insulating states with different occupation patterns, as well as supersolid states with 
different density distributions can be distinguished. 



I. INTRODUCTION 

Experiments with ultracold particles trapped in opti- 
cal lattices aim towards the engineering of exotic many- 
body quantum states Recently, the trapping and 
cooling of dipolar gases have attracted much attention 
0. The dipole moments induce long-range interactions 
between the particles, and new phases appear |3|. In 
the strongly correlated regime, it has been shown that 
there are many quasi degenerate metastable insulating 
states with denned occupation patterns 0-0|- These 
metastable states could be used for the storage and pro- 
cessing of quantum information in analogy to classical 
neural networks, where the information is robustly en- 
coded in the distributed stable states of a complex sys- 
tem [H, |j| . Another way to induce long-range interactions 
between atoms trapped in an optical lattice is via cou- 
pling to an external cavity mode. This has just recently 
been achieved experimentally and a checkerboard to a 
supersolid transition has been observed JlOJ . 

The detection of exotic strongly correlated phases re- 
quires novel experimental techniques that give access to 
high-order correlation functions. Proposals for detec- 
tion techniques typically make use of shot-noise mea- 
surements [Tl[ or atom-light interfaces [l2|]- Also, the 
counting statistics of atoms has been suggested as a tech- 
nique able to distinguish strongly correlated [U, [l4| and 
fermionic 11511 Hamiltonians, both at zero and finite tem- 
perature [16| . The detection of single atoms trapped 
in the o ptic al lattice has become experimentally avail- 
able jl7H20| only recently. Most counting experiments 
are performed after switching off the trapping poten- 
tial and letting the atoms propagate in the gravitational 
field. The counting statistics of Rb atoms falling within 
a high- finesse cavity has been reported in Ref. [21] . Also, 
fermionic and bosonic counting probability distributions 
have been measured for metastable Helium atoms falling 



onto a microchannel plate (22|, [23[ . 

The theoretical analysis of the counting process has 
so far mainly been considered for atoms trapped in the 
lattice. Propagation in the gravitational field mixes the 
initial modes of the atoms, such that the counting statis- 
tics in the lattice and after propagation are not expected 
to be the same. In this paper, we study the role of expan- 
sion in the counting process. We show that the mixing of 
the initial modes during the expansion becomes evident 
in the counting distribution when the detector is small 
compared to the size of the expanded wave function. We 
illustrate the effect by analyzing the counting statistics 
for bosons after time-of-flight expansion from the lattice. 
We consider initial states with different occupation pat- 
terns in the insulating regime and supersolid states with 
different density distributions in the superfluid regime. 
We calculate both the counting probabilities at a single 
detector and the joint probabilities at two detectors as 
a function of the horizontal distance between them. We 
show that a superfluid (SF) and Mott insulator (MI) state 
can be readily distinguished by their counting statistics. 
We further show that a suitable choice of the detector 
geometry allows for the detection of different occupation 
patterns in the insulating regime and different supersolid 
states. 

The paper is organized as follows. In Sec. [II] we re- 
view the propagation of the atomic wave functions and 
the atom counting formalism. In Sec. Mil we analyze the 
intensity of particles arriving at the detector, which con- 
sists of auto-correlation terms and crossed-correlations 
between the different expanded modes. Depending on 
the size and geometry of the detector, the ratio between 
the auto-correlations and the crossed-correlation terms 
changes. In Sec. IIVI we obtain closed expressions for the 
counting distributions for expanded superfluid and insu- 
lating bosonic states. We consider the counting statistics 
when using one detector and the joint counting distribu- 
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tion at two detectors. In Sec. [Vj we show our results 
and compare the SF with MI states and insulating states 
with different occupation patterns. 



II. DESCRIPTION OF THE SYSTEM 

We consider neutral bosonic atoms trapped in an opti- 
cal lattice. The system can be described using the Bose- 
Hubbard model [241 ] . which includes the hopping of the 
particles between neighbouring sites and the on-site two- 
body interactions. At zero temperature, the two limiting 
cases of the phase diagram are the SF state, where the 
hopping term dominates, and the MI state, where lo- 
cal interactions are dominant. The field operator 9(r,t) 
of the many-body system can be expanded into the N 
modes a* 



(1) 



For atoms trapped in an optical lattice, dj describes the 
destruction of a particle on site i. The corresponding 
initial wave functions are Wannier functions which are 
Gaussian functions centered at r, 



&(r,t = 0) 
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where the width u is chosen such that the initial wave 
functions at different sites i do not overlap. 

The atoms are released from the optical lattice and 
expand in the gravitational field. At finite t, we can 
apply the single-particle expansion 



0(r,t) 



j dr'K{r,r',t)<t>i{r\ 
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where the propagator for the free expansion in the grav- 
itational field reads 1251 



m \ 3 / 2 im(r 



2h 21m /i 
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The full propagated wave function is then written as 
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e 5RT^7 e s i^fT^) ; ( 5 ) 



where and r t = r + z t , with z t = (0,0, gt 2 /2) and we 
have used that \v t — r»| ^> to. Note that in the limit of 
u>t S> the expanded wave function is, up to a phase 
factor, a Gaussian function centered around z t with a 
width co t = ht/(mu)). 



A. Atom counting 

We describe a counting process in which the probabil- 
ity p{m) of counting to particles within a time interval r 



is measured at a detector located at a distance zq from 
the lattice. The probability of detecting to particles can 
be expressed as 0, [13] 



p{m) = 



(-iy 



-Q 



(6) 
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where the generating function Q(A) is given by the ex- 
pectation value of a normally ordered exponential of the 
intensity T, 

-xx 



C(A) = Tr(p 



(7) 



For photons, the intensity is proportional to an integral 
over the product of the negative-frequency part and the 
positive-frequency part of the field. The normal order- 
ing : ... : reflects the detection mechanism, in which the 
photons are absorbed at the detector, typically a photo 
multiplier or an avalanche photodiode. For the detection 
of atoms using microchannel plates, the detection process 
can be treated in an analogous way. 

Since typically not all the particles are counted, the 
intensity depends on the efficiency e of the detector and 
the detection time r. When the dynamics of the mea- 
surement are fast in comparison to the dynamics of 
the system, the intensity is proportional to the factor 
k = 1 — e _eT . For typical experimental situations, the 
dynamics of the system are determined by the expan- 
sion of the atomic cloud in the gravitational field, given 
by cot and the intensity can be described by the integral 
over the detector volume f2 of the positive-frequency and 
negative- frequency parts of the quantum fields describing 
the particles to be counted, multiplied by the efficiency 
factor k (H, 



1=K [ dr&(T,t d )*(T,t d ), 
Jq 



(8) 



where denotes the time at which the instantaneous 
measurement is performed. 

The formalism described above is easily generalized to 
the case of detection with multiple detectors [29| . For de- 
tection with M detectors, the generating function reads 

Q M (Ai, A 2 , .., A M ) = Tr(p : e~ £i A ^ :), (9) 

where the single detector intensity Xi for each of the de- 
tectors is given by eq. ([S]). For a configuration with two 
detectors, the joint probability distribution of counting 
to atoms at detector 1 and n atoms at detector 2 is given 
by 

ftfm+n d m+n 

P{m,n) = — — — Q 2 



d\fd)q^ Ai=l,A 2 =l' ^ 

We study the correlations corr(m, n) between the count- 
ing events detected at each detector by observing the 
ratio between the covariance and the single detector vari- 
ances, 

cov(to, n) 



coit(to, n) 



(11) 



a 2 (m)a 2 (n) ' 

where cov(m, n) = J2 m n mn P( m i n ) ~ "m* TO denotes 
the mean and a 2 (m) the variance of pirn). 
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III. DETECTION OF EXPANDING ATOMS 

Let us now discuss the counting process for the detec- 
tion of atoms expanding in the gravitational field. We 
consider a cubic detector located at a distance Zq from 
the lattice center with edge lengths A x , A y , A z . For sim- 



plicity, all through this paper we consider td = \J 2zq / g 
which is the time when the center of the cloud arrives 
at the detector. The intensity T of atoms registered at 
the detector defined in eq. © is thus determined by the 
expanded field operator of the atoms at the time td of 
detection, ^(zq, td)- Using eqs. (JT]) and (|8|), the intensity 
I takes the form 



where 



Aij(zo,n.,K) = n / dr(/)*(z a ,td)4'j(zo,td). 
Jn 



(12) 



(13) 



The elements of the correlation matrix Aij defined in eq. 
(fl~3|) describe the interference and autocorrelation terms 
between different modes registered at the detector. The 
diagonal terms represent the on-site correlations, whereas 
the off-diagonal terms represent the crossed-correlations 
between single particle modes initially located at different 
sites with distance \i — j\. 

Before studying the full counting distribution, let us 
consider the correlations given by the matrix elements 
A^. Using eq. (0 and assuming ui td ^> to, the autocor- 
relation elements are given by 



A; 



dr- 



1 



.3/2^3 



~2 — 



(14) 



For expanded wave functions at r > ri, the autocorre- 
lations become all equal and independent of the original 
lattice site i. The crossed-correlations are given by 



Aj, 



dr- 



1 



■3/2 



(15) 



The ratio between the crossed correlations and the auto- 
correlations depend crucially on the geometry of the de- 
tector. 

In Fig. [1] we show the on-site correlations eq. (| X4[) and 
the interference terms eq. (fT5)) in function of the size of 
the detector. We consider a one dimensional array in 
z-direction and plot the correlations at the location of 
the detector at (0,0, zq). We consider a fixed detector 
size in the x-y-plane, A = = A y and vary its width 
A z . Depending on the size of the detector, the whole 
cloud or a fraction of it is registered. For zq = 1 cm, 
the size of the expanded single-particle wave function at 
the detector is w td = 0.8 mm. In Fig. Q] a, we show 
that for detectors of size A z > 0.2 mm, the interference 
terms are negligible. This is easily understood from eq. 
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FIG. 1: The ratio between the on-site-correlations and inter- 
ference terms depend on the size of the detector, a) Wide de- 
tector limit A z ~ cot- b) Narrow detector limit A z <C Lot. We 
plot An (blue squares) Ai, i+ \ (green circles), A iti+ 2 (red di- 
amonds), and Ai,i+s (light blue triangles). Parameters used: 
zo — 1cm, A x = A y = 1cm and n = 1 



(Tj"4")) . as the auto-correlations are given by an integral 
over the detector volume around the center of a Gaussian 
function. For detectors that are large compared to the 
size of the cloud, the on-site correlations approach unity. 
In contrast, the interference terms eq. (|15p are given 
by the integral over a Gaussian function multiplied by 
a highly oscillating phase, such that they approach zero 
as the size of the detector increases. Thus, for on-site 
counting and for detectors which are larger than the size 
of the cloud, the crossed correlations disappear A^ ~ 
for i =/= j, while the auto-correlations approach An ~ k. 

As we show below, the detection of the auto- 
correlations between different modes is sufficient to dis- 
tinguish the long-range correlations in the system. In 
particular, we show that a MI state can be distinguished 
from a SF. 

On the contrary, as the auto-correlation terms for dif- 
ferent sites are equal, distinguishing states with different 
occupation patterns cannot be achieved in this limit. Fig. 
[T]b shows that for small detector sizes, the interference 
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terms are of the order of the on-site correlations. We will 
show that in this limit, different occupation patterns are 
distinguishable from the counting distribution. 

IV. ATOM COUNTING STATISTICS 

Let us now consider the counting distributions mea- 
sured at the detector after the expansion for different 
initial states of the system of atoms trapped in the lat- 
tice. 



A. Superfluid state 

First, let us focus on a SF state, ground state of the 
Bose-Hubbard model for very shallow lattices. We derive 
the counting distribution using the Gutzwiller ansatz (30l | 
for the wave function which assumes that it is a product 
of on-site coherent states. The initial state of the atoms 
in the lattice with N sites then reads: 

N 

m = n i«*>«' 

i 

where \(Xi)i is the coherent state on site i, 

oo n 

k), = e- M2/2 £^=N* (17) 

n— v 

and \n)i = (ai) ra |0) is a Fock state with n particles. Note 
that is an eigenstate of the annihilation operator 
^(r, t) of the expanded atoms, 

*(r,i)|V)=^^(r ) tK|V), (18) 

i 

where (pi is given by eq. ([2]). The state \ip) is thus 
an eigenstate of the expanded field operator ^(r,^) 
and we can write the generating function as Q(A) — 
e~^^ii a i a 3 A% 3 . Using eq. ((6]) the counting distribution 
p{m) reads 

(j2ij a i a j A ij) ^ , , 
p(m) = A i /_ e - E« (19) 

m! 

where Aij is given by eq. (fT"3|) . 

For a homogeneous superfluid with equal mean number 
of particles per sites, on — a for all i, and in the limit of 
big detectors where the diagonal elements of the matrix 
Aij are much bigger than the off-diagonal elements, the 
counting distribution of the SF simplifies to 

to! 

which corresponds to a Poissonian distribution with 
mean (and thus also variance) to = a 2 (m) = 



B. Mott Insulator state 

Let us now consider the Mott insulating regime. We 
first study a Mott insulator state with one particle per 
site, \ip) = 1 11.. 11). In this case, the generating function 
eq. ([7]) reads 

Q(X) = (ll..ll| : e-Wn*^ («•.**)*(«•.*-) : |H..H) 

= l - A J2 A u + x2 52( A u A jj + \ A a I 2 ) - (-2-1) 

i i<j 

We can rewrite eq. (|21[) using the minors of the matrix 

N 

Q(A) = l + ]T(-l) fc A fc M + (Afc), (22) 

k=l 

where M + (^4,to) denotes the permanent perm(A) = 
So-eS ^i=i A i,cr(i) of the corner blocks of size to of the 
matrix A. Note that M + (A,k) is closely related to the 
principal minors of the matrix, which are defined as 
the determinant of the respective block matrices. The 
counting distribution p(m) can then be calculated using 
Eqs. © and |22j) . 

As was outlined above, in typical experimental situa- 
tions the detector is far away from the lattice and much 
bigger than the cloud, such that the off-diagonal elements 
of A^ are negligible and the diagonal elements An are 
equal for all i. In this case the generating function Q for 
the Mott insulator state with unit filling is given by 

Q ( A ) = E (J (~ A ^) fe = i 1 - A ^ N > ( 23 ) 

where Ad denotes any of the (equal) diagonal elements. 
The counting distribution p(m) is then given by 

p(m)= ( N )Af(l-A d ) N - m (24) 

\TO / 

This corresponds to the distribution of a fock state. The 
mean m and variance a 2 (to) of the distribution are given 
by 

m = NA d , a 2 (m) = NA d {l-A d ) (25) 

Let us now consider the different occupation patterns 
that arise in the strongly correlated regime. In particu- 
lar, we focus on such states where at most one particle 
occupies each site. The generating function is then calcu- 
lated by eq. (f2"2")l . with a correlation matrix A', composed 
of the elements of the correlation matrix A in eq. (fT5|) 
multiplied by the occupation numbers n-i and rij of the 
involved sites, 

A' a, ii. A,,. (26) 

Finally, let us consider a symmetric superposition of 
all possible states with filling factor N p /N s , where N p is 
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the number of particles, N s denotes the number of sites 
and N p < N s , the generating function reads 

2 = 1 + Y,(-l) m \ m F M i(A, m, N p: N s ), (27) 

m 

where 

(N B -m\ 

F MI (A,m,N p ,N s ) = S^^L M +(A,m) 

\Np) 

+ -^--2" l /C(m), (28) 

\N P ) 

where /C(m) is defined as the mfold product over the 
sum with non- repeated indices of the real part of Aij, 
EiKjMAii) ■ For m = 2, e.g. M+ = E,<,(^% + 
\Aij\ 2 ) and K{m) = Re(Aij)Re(Aki) with k,l ^ 



C. Counting at two detectors 

In this section, we consider the detection of the MI and 
SF state using two detectors and study the correlations 
between the counting events. 

For the MI state, the joint counting distribution 
p(m, n) of counting m particles at one detector and n 
particles at the other is given by eq. (|10[) . where the 
generating function for two detectors is given by 

N 

Q 2 = ^(-l) fe M+(A 1 A( 1 ) + X 2 A^,k). (29) 

k=l 

For detectors that are located symmetrically with respect 
to the origin in the ir-y-plane, in typical experimental 
situations the off-diagonal elements of A^ are negligible 
(see fig. QJ, and the diagonal elements Ad are all equal 
for both detectors, = A d 2 ' = Ad- The generating 
function thus simplifies to 

Q 2 (A 1 ,A 2 ) = £r)(- J 4 rf ) fe (A 1 + A 2 ) fc 

fc=0 ^ ' 

= (l-Ai(A 1+ A 2 )) w , (30) 

and the counting distribution is given by 

p{m,n) = (-l)"+ m (l - 2A d ) N P- m - n x 

{-A d ) m+n , ,," p! Yf (31) 

For the SF state, the joint counting distribution 
Psf(w, n) is the product of the two single detector distri- 
butions pi (m) and p 2 in) given by eq. (1191) . The counting 
events at the two detectors are thus not correlated. 
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FIG. 2: MI vs. SF as a function of distance from detector. 
Probability distribution for MI (black bars) and superfluid 
(white bars) states in a 3x3x3 lattice. A x = A s = 2 mm, 
A z = 2 cm, k = 1. a) zq = 1 cm, b) zq — 3 cm, c) zq = 5 cm 



V. RESULTS 

A. Mott Insulator and Superfluid state 

We consider the counting distributions of a SF and a 
MI state of bosons with the same average number of par- 
ticles released from a three dimensional optical lattice. 
We assume the limit of a large detector, where the count- 
ing distribution is determined by the on-site correlation 
terms. In Fig. [5J we plot the counting distributions for a 
SF and a MI state at different distances between the de- 
tector and the lattice. With increasing distance from the 
detector, a smaller fraction of the expanded wave func- 
tion is registered. The difference between the MI and the 
SF becomes less visible, and the mean of the counting 
distribution decreases. In Fig. [31 we plot the mean and 
the variance of the counting distributions, both normal- 
ized by dividing by N, for the superfluid and the Mott 
insulator state for a detector with fixed size at different 
distances z$ from the lattice. 



B. Mott Insulator and Superfluid state with two 
detectors 

Let us now consider two detectors of the same size that 
are placed symmetrically at a distance xi = (x<j,0, Zo) 
and x 2 = (— Xd-, 0, zq) from the lattice center. In the 
limit of large detectors, we study the joint counting dis- 
tribution of the SF and MI state for different distances 
between the detectors. Fig. [4] shows the counting distri- 
butions for two overlapping detectors (left column) and 
for two detectors separated by 2x d = 1cm (right col- 
umn). For the SF state, shown in the lower row in figs. 



Z Q (cm) 

FIG. 3: a 2 (m)/N of the counting distribution for the MI 
(blue squares) and SF (green circles) state with respect to 
the distance from the detector Zo. A x = A y = 2 mm, A z — 2 
cm, k = 1. 

IH the joint counting distribution is a Gaussian function 
for both cases. This is expected, as the joint counting 
distribution eq. (|T5)) is a product of the single detector 
counting distributions. This is analogous to the detec- 
tion of coherent states of light. For the MI state shown 
in the upper row in fig. |4j we observe a squeezed distri- 
bution, indicating the correlations of the atoms counted 
at the two detectors. Note that as the distance between 
the detectors increases, the squeezing of the distribution 
is less pronounced. The correlations between the count- 
ing events at the two detectors can be seen more clearly 
when looking at the correlation function eq. (fTTj) . Note 
that for the superfluid state, there is no difference be- 
tween the joint counting distribution and the product of 
the single particle distributions. For the Mott state, we 
study the correlations for varying distance between the 
two detectors Xd- In Fig. [5j we show how the correlations 
decrease when increasing the distance between detectors 
Xd- Note that the distance Xd denotes the distance be- 
tween the center of the two detectors. For Xd = 0, the 
detectors fully overlap, and for Xd > A the detectors are 
completely separated. 

C. Detection of insulating states with different 
occupation patterns 

Let us now focus on the detection of insulating states 
with different occupation patterns by particle counting. 
As discussed above, in order to detect the different pat- 
terns, the crossed-correlations have to be of the order of 
the autocorrelations. This is clear as away from the lat- 
tice, all the on-site correlation terms become equal. Let 
us discuss the example of a checkerboard state, where 
every second site is occupied, and a state with stripes, 
where every second line is occupied. For the striped state, 
the leading crossed-correlation terms eq. (fT5)l are the 
ones that correspond to the nearest neighbors. For the 
checkerboard state, where neighboring sites are not oc- 
cupied, the leading terms are the ones that correspond to 
diagonally adjacent sites. In order to distinguish the dif- 
ferent patterns, it is thus essential that these two leading 
crossed-terms are sufficiently different and at the same 



m o m 



c > d) 




FIG. 4: Joint probability distribution of an expanded MI (up- 
per row) and a SF (lower row) in a 4x4 lattice with two sym- 
metrically placed detectors. In fig. a) and c) xa = 0. In fig. 
b) and d) Xd = 1 cm. Parameters used: zo = 1cm, A z = 2 
mm, A x — A y — 2 cm, k — 0.5. 




FIG. 5: Correlations of the joint probability distribution for 
an expanded MI state with two symmetrically placed detec- 
tors. As the distance between the detectors increases, the 
counting events at the two detectors are no longer correlated. 
Parameters used: zo = 1 cm, A z = 2 mm, A x = A y = 2 cm, 
k = 0.5 



time comparable to the on-site correlations. From Fig. 
[TJ we see that this implies that the limit of small detec- 
tors has to be considered. However, if the detector is 
very small, all the terms are equal and the patterns are 
not distinguishable. One should thus consider an inter- 
mediate detector size. 

In Fig. O we illustrate the effect for a ID system of 
N = 12 particles. We compare the counting distribu- 
tions of a checkerboard-like state, where every second 
site is occupied, and a state where a block of six sites 
is occupied and a block of six sites is empty. In order 
to distinguish the two states, from Fig. [TJ we choose a 
detector size of A = 0.02 mm, such that the ratio of the 
crossed-correlation terms between neighboring sites and 
the autocorrelations is 0.6. Fig. [B] shows that the dif- 
ferent occupation patterns are reflected in the counting 
distribution. 
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FIG. 6: The counting distributions of an expanded one- 
dimensional checkerboard (black bars) and striped insulating 
pattern (white bars) are clearly distinguishable. Parameters 
used: zq = 1 cm, A = 0.1 cm, A z = 0.02 mm,K = 1, N p = 12 



D. Detection of a Supersolid state 

As for the detection of states with different occupa- 
tion patterns in the insulating regime, the detection of 
supersolid states [U HH requires the limit where the 
crossed-correlation terms for neighboring sites are com- 
parable to the auto-correlation terms. We consider a su- 
persolid state with N sites and mean density ct^i = and 
a%i-\ = 7. For the limit where the crossed-correlation 
terms for neighboring sites are the only non-negligible 
interference terms, the counting distribution eq. (1191) is 
given by a Poissonian distribution with mean 

m= ^-A d {(3 2 + 1 2 ) + 2NA NN f3 1 , (32) 

where A d denotes the diagonal elements corresponding 
to the on-site correlations and Ann denotes the nearest 
neighbor crossed-correlation terms. Let us compare this 
to a superfluid state with a homogeneous density per site, 
l^.j 2 _ \@\ +n f or a u ^ rp^g coun ^j n g distribution eq. 
(I19p is thus given by a Poissonian distribution with mean 

m = ^A d ((3 2 + 7 2 )+NA NN (0 2 +~/ 2 ). (33) 

From eqs. (|32j) and (f33|) it is clear that a supersolid state 
can be distinguished from a superfluid state by particle 
counting. In Fig. [7] we illustrate this by comparing a 
supersolid state to a superfluid state. 

VI. SUMMARY 

We have studied the counting distributions of atoms 
falling from an optical lattice and propagating in the 
gravitational field. The intensity of atoms recorded 
at a detector located far away from an optical lattice 



can be decomposed into autocorrelation and crossed- 
correlations between the expanding modes. The ratio 
between these terms depends crucially on the geometry 
of the detector. In the limit when the detector is large 
compared to the expanded modes, the crossed-correlation 
terms are negligible and only long-range correlations of 
0.1 1 . 1 1 . 1 1 




m 



FIG. 7: The counting distributions of an expanded supersolid 
state with \/3\ 2 — 0.5 and \j\' 2 = 1.5 (black bars) and a super- 
fluid state |a| 2 = 1 (white bars) are clearly distinguishable. 
Parameters used: zo = 1 cm, A = 1 cm, A z = 0.02 mm, 
K = 1 



different states can be distinguished. In this limit a SF 
state has a poissonian number distribution while a MI 
has subpoissonian number distribution for a detector of 
finite size located at a distance zq from the lattice. The 
two states can also be readily distinguished from the joint 
probability distribution of counting the particles at two 
detectors. In the SF regime, the joint probability dis- 
tribution is a product of the two independent number 
distributions while in the MI regime, the distributions 
are highly correlated. 

When the detector is small compared to the expanded 
wave function, the crossed-correlation terms for adjacent 
sites are of the order of the auto-correlations. We have 
shown that by choosing the size of the detector in an ap- 
propriate way, different occupation patterns can be dis- 
tinguished by particle counting after expansion both in 
the insulating as well as in the superfluid regime. 
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